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Abstract 

Next-to-leading order calculation for three jet heavy quark production in e + e~ collisions, 
including complete quark mass effects, is reviewed. Its applications at LEP/SLC are also 
discussed. 
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Next-to-leading order calculation for three-jet heavy quark production in e + e~-collisions, including complete 
quark mass effects, is reviewed. Its applications at LEP/SLC are also discussed. 



The importance of the corrections due to the 
mass of the heavy quark in the jet-production in 
e + e _ -collisions has been already seen in the early 
tests of the flavour independence of the strong 
coupling constant The final high precision 

of the LEP/SLC experiments required accurate 
account for the bottom-quark mass in the the- 
oretical predictions. If quark mass effects are 
neglected, the ratio a b s /a^ ds measured from the 
analysis of different three-jet event-shape observ- 
ables is shifted away from unity up to 8% || (see 
also |). 

Sensitivity of the three-jet observables to the 
value of the heavy quark mass allowed to con- 
sider the possibility of the determination of 
the b-quark mass from LEP data, assuming the 
universality of a s . In a recent analysis of three- 
jet events, DELPHI measured the mass of the 
b-quark, mj, for the first time far above the pro- 
duction threshold Q. This result is in a good 
agreement with low energy determinations of nib 
using QCD sum rules and lattice QCD from T 
and £?-mesons spectra (for recent results see e.g. 
|7|) The agreement between high and low energy 
determinations of the quark mass is rather im- 
pressive as non-perturbative parts are very dif- 
ferent in the two cases. 

In this contribution we will discuss some as- 
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pects of the next-to-leading order (NLO) calcula- 
tion of the decay Z — > Sjets with massive quarks, 
necessary for the measurements of the bottom- 
quark mass at the Z-peak. Recently such cal- 
culations were performed independently by three 
groups §|(j . 

The first question we would like to answer 
whether it is not at all surprising that LEP/SLC 
observables are sensitive to mi, as the main scale 
involved is the mass of the Z-boson, Mz 3> 
Indeed, the quark-mass effects for an inclusive 
observable such as the total width Z — > bb are 
negligible. Due to Kinoshita-Lee-Nauenberg the- 
orem such observable does not contain mass sin- 
gularities and a quark-mass appears in the ratio 
m 2 q (Mz)/M z w 10~ 3 , where using MS running 
mass at the M^-scale takes into account the bulk 
of the NLO QCD corrections Jll|,§. 

However, the situation with more exclusive ob- 
servables is different. Let's consider the simplest 
process, Z — > bbg, which contributes to three-jet 
final state at the leading order (LO). When the 
energy of the radiated gluon approaches zero, the 
process has an infrared (IR) divergence and in 
order to have an IR-finite prediction, some kine- 
matical restriction should be introduced in the 
phase-space integration to cut out the trouble- 
some region. In e + e~-annihilation that is usually 
done by applying the so-called jet-clustering algo- 
rithm with a jet-resolution parameter, y c (see Jl2) 
for recent discussion of jet-algorithms in e + e~). 
Then the transition probability in the three-jet 



part of the phase-space will have contributions as 
large as l/y c ■ (ml/M z ), where y c can be rather 
small, in the range 10~ 2 — 10~ 3 . Then one can ex- 
pect a significant enhancement of the quark-mass 
effects, which can reach several percents. 

The convenient observable for studies of the 
mass effects in the three-jet final state, proposed 
some time ago 0H, is defined as follows 



R b 3 d = 



r| J -( 2 / c )/r 6 



r d 3j (yc)/T d 
l + r& (bo{yc,n) 



(l) 



where F^ and T q are three-jet and total decay 
widths of the Z-boson into quark pair of flavour 
q, r& = m|/M§. Note that above expression is 
not an expansion in r&. 

The LO function, bo, is plotted in fig. f for four 
different jet algorithms. 
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Figure 1. LO contribution to the ratio R 
function of y c ( see eq. for the definition ) for 
nib = 3GeV (dashed curve) and m& = 5GeV 
(solid curve). 

Together with well known JADE, E and 
DURHAM schemes we consider the so-called EM 
algorithm || with a resolution parameter j/y = 
2piPj/s and which was used for analytical calcu- 



lations in the massive case [[| . The main observa- 
tion from fig. 1 is that for y c > 0.05, bo is almost 
independent of the value of mb for all schemes. 
Although, this remains true also for smaller y c in 
DURHAM and E schemes, there is a noticeable 
mass dependence in JADE and EM schemes. 

Note that bo is positive for E-scheme. That 
contradicts the intuitive expectations that a 
heavy quark should radiate less than a light one. 
This unusual behavior is due to the definition 
of the resolution parameter in E-scheme, y^ = 
(j>i +Pj) 2 1 s, which has significantly different val- 
ues for partons with the same momenta in the 
massive and massless cases, and it can be used as 
a consistency check of the data. 

In what follows we restrict ourselves to 
DURHAM scheme, the one used in the experi- 
mental analysis B], and bo can be interpolated 



as: bo 



Jo) 



b { o 2) ln 2 y c 



In the LO 



calculations we can not specify what value of the 
b-quark mass should be taken in the calculations: 
all quark-mass definitions are equivalent (the dif- 
ference is due to the higher orders in a s ). One 
can use, for example, the pole mass Mb ~ 4.6GeV 
or the MS'-running mass rh~b(p) at any scale rel- 
evant to the problem, mb < fx < Mz, with 
m b (m b ) « 4.I3GeU and M b (M z ) w 2.83GeU. As 
a result, the spread in LO predictions for differ- 
ent values of b-quark mass is significant, the LO 
prediction is not accurate enough and the NLO 
calculation should be done. 

At the NLO there are two different contribu- 
tions: from one-loop corrections to the three- 
parton decay, Z — > bbg and tree-level four-parton 
decay, Z — > bbgg and Z — > bbqq, q = u,d, s, c, b 
integrated over the three-jet region of the four- 
parton phase-space. In the NLO calculation one 
has to deal with divergences, both ultraviolet 
and infrared which appear at the intermediate 
stages. The sum of the one-loop and tree-level 
contributions is, however, IR finite. We would 
like to stress that the structure of the NLO cor- 
rections in the massive case is completely different 
from the ones in the massless case lfl3| . That is 
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due to the fact that in the massive case, part of 
the collinear divergences, those associated with 
the gluon radiation from the quarks, are softened 
into lnrt and only collinear divergences associ- 
ated with gluon-gluon splitting remain. 

In the NLO calculations one should specify 
the quark mass definition. It turned out tech- 
nically simpler to use a mixed renormalization 
scheme which uses on-shell definition for the 
quark mass and M S definition for the strong cou- 
pling. Therefore, physical quantities are origi- 
nally expressed in terms of the pole mass. It can 
be perfectly used in perturbation theory, however, 
in contrast to the pole mass in QED, the quark 
pole mass is not a physical parameter. The non- 
perturbative corrections to the quark self-energy 
bring an ambiguity of order rj 300MeV (hadron 
size) to the physical position of the pole of the 
quark propagator. Above the quark production 
threshold, it is natural to use the running mass 
definition (we use MS). The advantage of this 
definition is that m&(/u) can be used for /i 3> mj. 
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Figure 2. NLO function b\ for different mb (see 
eqs.^),^) for the definition and text for details). 
The errors are due to numerical integrations. 

The pole, Mb, and the running masses of the 
quark are perturbatively related 



M b = mb(fj) 



We use this one-loop relation to pass from the 
pole mass to the running one, which is consistent 
with our NLO calculations. To match needed pre- 
cision we have to use this equation for values of /i 



about nib- Then we can use one- loop renormaliza- 
tion group improved equation in order to define 
the quark mass at the higher scales. Substituting 
eq.(H) into definition eq.(Jj]) we have 



R b 3 d (y c ,m b (n),ti) = 



l + r b {p) 



a s (fj,) 



2 & o m^f 

A* 



(3) 



with bi = b x +& (8/3-21nr 6 ) and r b =ml/M%. 

In fig. 2 we show the NLO function bi(yc,n) 
calculated for three different values of the quark 
mass : 3GeV (open circles), 4GeV (squares) and 
5GeV (triangles). 
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Figure 3. The ratio R\ d (eq. ^)). Solid curves 
- LO predictions, dashed curves give the NLO re- 
sults (see text for details). 

In contrast to the bo, one sees a significant 
residual mass dependence in b\, which can not 
be neglected. The solid lines in fig. 2 represent a 
fit by the function: b\ = b^ + b^ hiyc + bf ^ lnrf, 
performed in the range 0.01 < y c < 0.1 The qual- 
ity of this interpolation is very good and the main 
residual nib dependence in b\ is taken into account 
by In rb term. Inclusion of higher powers of In rb 
does not improve the fit. 

Fig. 3 presents theoretical predictions in the 
DURHAM scheme for the R\ d observable mea- 
sured by DELPHI 0. The solid lines are 
LO predictions for the b-quark mass, mb = 
rh~b{Mz) — 2.83GeV (upper curve) and mb = 
Mb = 4:.6GeV (lower curve). The dashed curves 
give NLO results for different values of scale 
(i : 10, 30, 9lGeV. One sees that NLO curve 



for large scale is naturally closer to LO curve for 
mb(Mz) and for smaller scale is closer to the LO 
one with m& = Mi,. 

Fig. 4 illustrates the scale dependence of R^ d 
for y c — 0.02. By studying the scale depen- 
dence, which is a reflection of the fixed order 
calculation, we can estimate the uncertainty of 
the predictions. The dashed-dotted curve gives 
//-dependence when eq. (Q) was used, so it is /i- 
dependence due to renormalization of the strong 
coupling constant, a s . 
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Figure 4. The ratio R^ d as a function of the scale 
M for y c = 0.02. 

Other curves show //-dependence when R^ d 
is parameterized in terms of the running mass, 
mb(Mz), eq.(|^), but different mass definitions 
have been used in the logarithms. The conserva- 
tive estimate of the theoretical error for the R^ d 
is to take the whole spread given by the curves. 
The uncertainty in R\ d induces an error in the 
measured mass of the b-quark, AR^ d = 0.004 — * 
Ami, — 0-23GeV. This theoretical uncertainty 
is, however, below current experimental errors, 
which are dominated by fragmentation. 

To conclude, the NLO calculation is necessary 
for accurate description of the three-jet final state 
with massive quarks in e + e _ -annihilation. Fur- 
ther studies of different observables and different 
jet-algorithms could be very useful for the reduc- 
tion the uncertainty of such calculation. 
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